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I present a study of the equation of state in quenched QCD, discussing some systematic effects re- 
lated to the lattice geometry. In particular, I comment on the modification of the Stefan-Boltzmann 
law for a gas of free gluons in a finite system, and study the impact it might have on numerical 
results at high temperatures, for the typical parameters of current lattice simulations. Finally, I 
apply the results of this study to the analysis of data obtained from simulations of SU(N) gauge 
theories with N > 3 colors, in a temperature range up to 3T C , where infrared effects appear to be 
under control. Preliminary results for various thermodynamic observables for SU(4), SU(5) and 
SU(6) gauge theories are found to be close to each other and to those for SU(3), in agreement 
with other similar studies. This may be relevant for the theoretical description of the QCD plasma. 
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1. Introduction and motivation 



Experiments at temperatures J up to a few times the deconfinement temperature T c reveal a 
strongly interacting system [|l|], which behaves as a nearly perfect fluid [0]. In this regime, lattice 
calculations are essentially the only tool to determine the equation of state (EoS) from the first 
principles of QCD. 

In particular, for the gluon sector of the QCD plasma, the lattice results show strong deviations 
from the Stefan-Boltzmann (SB) limit, with a large deficit in the entropy and pressure pi]. This led 
to conjecture [Q] that the gluon plasma may admit an effective description based on the AdS/CFT 
correspondence [Q]. Together with recent efforts to apply similar gauge/gravity techniques to build 
a holographic dual of QCD this has triggered interest in lattice studies of SU(Af) Yang-Mills 
thermodynamics wifh./V > 3 colors [0-^]. 

One aspect of the SU(iV > 3) simulations is that, in general, they can be carried out using 
smaller lattices than those needed for SU(3) — see, e.g. refs. [|Io|]. This should not pose a problem 
of infrared (IR) effects, because for temperatures of the order of a few times T c finite-volume effects 



are expected to be exponentially suppressed, due to the existence of screening masses [11]. 



However, it was recently suggested [ |12| ] that nontrivial IR corrections affecting a gas of free 
gluons may still be relevant for the strongly interacting plasma at temperatures 0{T C ). For a pe- 
riodic box of volume L and timelike size l/T, these corrections depend logarithmically on the 
aspect ratio of a timelike cross-section x := LT, and vanish for x — > oo. In this contribution I study 
numerically the impact of these effects at temperatures up to about 3T C . Issues related to lattice IR 



effects were also studied in refs. [13]. 



2. Finite- volume corrections to the partition function of the free gluon gas 

In the continuum, the partition function ^ for a gas of free SU(jV) gluons in a finite box of 
spacelike sizes L x L x L at temperature T was calculated exactly in ref . []12[] : 
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The leading finite- volume correction to lni2° is a logarithmic function of LT. Neglecting terms 
0(e kLT ), the pressure p, energy density e, and free energy density / in a finite volume read: 
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Note that p is no longer equal to — /, which is in contrast with the usual assumption underlying 



the determination of the EoS on the lattice with the integral method []14|]. The accuracy of these 
corrections for the T — > °° limit on the lattice is manifest when comparing them with, e.g., the 
finite- volume corrections to e obtained by numerical integration in ref. [JT3J] — see fig. 1 in ref. [|l2|]. 

One may wonder whether these effects can also play a role for a gas of strongly interacting 
gluons at finite temperature — in particular, at relatively low temperatures 0(T C ). As the logarith- 
mic corrections appearing on the r.h.s. of eq. ( |Q| ) stem from the regularization of the divergent 
contribution to the functional integral coming from constant configurations, and periodic b.c. allow 
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the existence of constant configurations also in the interacting system, it may be that a nontrivial IR 
correction still affects the gluon gas at relatively low temperatures. In particular, in ref. [12] it was 
pointed out that, for the typical parameters of present lattice simulations, this effect could account 
for a large fraction of the deviations from the SB limit observed at temperatures of the order of T c . 



3. Numerical results 

I run simulations of S\J{N) gauge theories with the standard isotropic Wilson action, using an 
algorithm that combines heat-bath updates [ |I6| ] for SU(2) subgroups JT7| ] and full-SU(A r ) overre- 



laxation updates [18]; part of the simulations at T = were run using Chroma [|19|]. The lattice 
extent in the timelike direction in lattice units was N t = 5 (and, in some cases, 6); the spacelike 
volume in lattice units was N 3 , with N s up to 22, 18, 16 and 16 for SU(3), SU(4), SU(5) and 
SU(6), respectively. The T = simulations were run on N* lattices. 

For SU(3), I set the scale using r [20], while for the other groups I used the string tension [pl|]. 



I determined the EoS using the integral method [14], from measurements of the average plaquette 



at finely separated j8 values: in particular, n = 180 intervals of amplitude A/3 = 0.005 were used 
for SU(3). The numerical integration was done with the method described by eq. (A.4) in ref. [p2|], 
where it was used to measure the interface tension in the Ising model 1 ; the errors of this method 
are 0(n~ 4 ). Cutoff effects on the asymptotic normalization of the EoS were corrected using the 
Ri(N T ) factor g^]. 

The IR corrections to the EoS can be calculated analytically only for the free-gluon gas; the 
result obtained in this limit can be considered as an upper bound for possible nontrivial IR effects 
in the strongly interacting theory at temperatures of the order of T c . To check if in this regime the 
system is sensitive to logarithmic finite-volume corrections, I calculated the pressure from lattices 
characterized by different values of LT in two different ways. First, assuming that the data are not 
affected by nontrivial IR effects. Second, assuming maximal sensitiveness, i.e., assuming the same 
corrections as for the free-gluon gas. 2 If the IR corrections affecting the gas of free gluons were 
still relevant at T ~ 2T C or 3T C , then the latter method should give better consistency among results 
obtained from different lattices. 

This, however, seems not to be the case: the left panel of fig. [I] shows that the discrepancies 
between the results obtained using the first method are compatible with statistical fluctuations; this 
still holds up to T ~ 3T C (central panel). On the contrary, the second method overcompensates 
the differences between the two data sets, driving the curves away from each other (right panel). 
Similar results also hold for SU(4). 

Therefore, within the precision of our data, at temperatures up to T ~ 3T C we do not observe 
evidence for the logarithmic IR corrections that affect the free-gluon gas. Yet, as these corrections 
should eventually show up in the perturbative regime at high enough temperatures, it would be in- 
teresting to study at what T values they can be seen on the lattice. Since the logarithmic corrections 
are nonnegligible only for LT = 0(1), this problem could be studied by running lattice simulations 



1 See refs. [ p3| ] for further details. 

2 Since the quantity that is measured on the lattice is — /, rather than p, the latter method requires that the data be 
shifted (to compensate for p + f ^ 0) and rescaled (to account for the different high-temperature limit). 
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Figure 1: Left panel: The EoS for SU(3), evaluated neglecting possible logarithmic corrections to the 
partition function, for lattices of size 20 3 x 5 (black crosses) and 18 3 x 5 (red circles). The plot shows 
the results for p/T 4 (normalized to the SB limit) vs. T JT C . Central panel: High-temperature zoom of 
the previous plot. Right panel: The EoS evaluated assuming maximum sensitiveness to IR effects (i.e. 



compensating for the logarithmic effects that affect a free-gluon gas), according to eq. (2.2). 




Figure 2: The rescaled pressure p/T 4 (normalized to the SB limit) from simulations on a 16 3 x 5 lattice, vs. 
T/T c , for SU(4) (left panel), SU(5) (central panel) and SU(6) (right panel). 



at large j8 values and shrinking physical volume. 3 Although such a limit has its own reasons of 



theoretical interest [27], in the present work I restricted my analysis to temperatures 0(T C ), and to 
sufficiently large physical volumes. 

Having found no evidence for logarithmic IR corrections in the T <3T C regime, I proceeded 
to the evaluation of thermodynamic observables in SU(iV). Figs. |2||5] show preliminary results for 
equilibrium thermodynamics observables (pressure, trace of the stress tensor A = e — 3p, energy 
density and entropy density s), for 3 < N < 6, normalized to their SB limits. In agreement with the 
conclusions of ref. [^], the results from various gauge groups are essentially compatible with each 
other (with small quantitative corrections), and clearly different from the SB limit. Further results, 
including data for the SU(8) gauge group, are presented elsewhere [^J. 



3 The SU(3) EoS at very high temperatures was studied in ref. [p5[, where the results shown in fig. 5 were obtained 
from simulations with N t = S and A' 3 = 24 3 or iV 3 = 32 3 [ ^o] | . The maximal IR correction would then shift the asymptotic 
value down to about 97% of the SB limit, which indeed looks compatible with their result at the highest temperature 
T = 3 x 10 7 r e (although the data are also compatible with the conventional SB limit). Note however, that this comparison 
is only approximate, because in the numerical method of ref. J25| ] the geometric parameter x is not held fixed, and 
furthermore the pressure at high temperatures is obtained using the integral method, starting from low temperatures. 
This is based on the p = —f equality, which — as we have seen — is reliable in the low-T regime, but would also start 
being affected by IR corrections when T is increased. 
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Figure 3: Left to right: A/T 4 , normalized to £{N 2 - 1), vs. T/T c , for SU(3), SU(4), SU(5) and SU(6). 




Figure 4: Left to right: The rescaled energy density e/T 4 (normalized to its SB limit) vs. T/T c , for SU(3), 
SU(4), SU(5) and SU(6). 




Figure 5: Same as in fig. Q but for the rescaled entropy density s/T 3 



4. Summary and outlook 



The nontrivial IR corrections affecting a gas of free gluons in a finite box were derived in 



ref. [|12p, where it was also suggested that they may still play a role at temperatures of the same or- 
der of magnitude as T c . To the level of precision achieved, my numerical results in a (phenomeno- 
logically relevant) temperature range up to about 3T C did not reveal any nontrivial finite-volume 
effects. However, the latter might probably be observed in numerical simulations at very high tem- 
peratures and small volumes [|25|]. Next I also measured bulk thermodynamic quantities in S\J(N) 
gauge theories with 3 < N < 6 colors, at temperatures up to 3T C . In agreement with ref. 0, the 
results for different S\J(N) groups share similar features, and the strong deviations from the SB 
limit observed in SU(3) survive the large-Af limit. This leans support to the hypothesis that QCD 
may admit a nearly conformal effective description at temperatures of a few times T c , where the 



theory is still strongly interacting [28]. 
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